Curve crossing induced dissociation : An analytically solvable model 

Aniruddha Chakraborty 

School of Basic Sciences, 

Indian Institute of Technology Mandi, 

Mandi, Himachal Pradesh, 750001, India. 

(Dated: December 26, 2011) 

In our earlier papers we have proposed an analytically solvable model for the two state curve crossing 
problem which assumes the coupling to be a Dirac delta function. It is used to calculate the effect 
of curve crossing on electronic absorption spectrum and Resonance Raman excitation profile for 
the case of harmonic potentials. In this paper we have extended our model to deal with the curve 
crossing induced dissociation cases. Our method is used in this paper to calculate the effect of curve 
crossing induced dissociation on electronic absorption spectrum and Resonance Raman excitation 
profile. In this paper, a model consisting of a Harmonic oscillator and a Morse oscillator, coupled 
j.-^ ' by Dirac delta function, is solved. 
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INTRODUCTION 



■"nJ" i Nonadiabatic transition due to potential curve crossing is one of the most important mechanisms to effectively induce 

" ■ electronic transitions in collisions |ll4l3l|. This is a very interdisciplinary concept and appears in various fields of 

'_ physics, chemistry and biology [2|. Two state curve crossing is in general classified into the following two cases 

,^ ' according to the crossing scheme: (1) Landau-Zener (L.Z) case, in which the two diabatic potential curves have the 

^-i[ same signs for the slopes and (2) non-adiabatic tunnelling (N.T) case, in which the diabatic curves have the opposite 

sign for slopes. There is also a non-crossing non-adiabatic transition called the Rosen-Zener-Demkov type [l|, y]) in 



which two adiabatic potentials are in near resonance at large R. The theory of non-adiabatic transitions dates back 



^ to 1932, when the pioneering works for curve-crossing and non-crossing were published by Landau J14i] . Zener [IE 



Q^[ and Stueckelberg [16| and by Rosen and Zener [17| respectively. Since then numerous papers by many authors have 
been devoted to these subjects, especially to curve crossing problems [l|, l^]. In our earlier papers we have proposed 
an exactly solvable model for the two state curve crossing problem which assumes the coupling to be a Dirac delta 



J> ■ function |8l-lll| . This model is used to calculate the effect of curve crossiiig on electronic absorption spectrum and on 



CN ' Resonance Raman excitation profile for the case of harmonic potentials |8l-lllj . We have later generalized our model 
\l to deal with general multi-channel curve crossing problem too ,12;]. Even very recently our model ia extended to deal 

^-j I with nonadiabatic tunneling in an ideal one dimensional semi- infinite periodic potential swtems [13|. We have also 

proposed an analytical method for the two state curve crossing problem for any coupling llq]. We have also used our 
t^ \ analytically solvable to deal with scattering problems 19']. The same method has been applied recently to the case 

of predissociation [20J . Our work is in progress to deal with nonadiabatic tunneling in an ideal one dimensional finite 



periodic potential systems 2l|. In this paper we have extended our model to deal with the curve crossing induced 
dissociation cases. In this paper we analyze the effect of curve crossing induced dissociation on electronic absorption 
spectrum and Resonance Raman excitation profile using a model consisting of a Harmonic oscillator and a Morse 
oscillator, coupled by Dirac delta function. 

THE MODEL 

We consider two diabatic curves, crossing each other. There is a coupling between the two curves, which causes 
transitions from one curve to another. This transition would occur in the vicinity of the crossing point. In particular, 
it will occur in a narrow range of cc, given by 

\V^{x)-V2{x)\:^\V{x,)\. (1) 

where x denotes the nuclear coordinate and Xc is the crossing point. Vi and V2 are the diabatic potentials and V 
represent the coupling between them. In reality the transition between Vi (x) and V2 {x) occur most effectively at the 
crossing, because the necessary energy transfer between the electronic and nuclear degrees of freedom is minimum 
there. Therefore it is interesting to analyze a model, where coupling is localized in space near Xc rather than using a 
model where coupling is same everywhere (i.e. constant coupling). Thus we put 

V{x)=Ko5{x-x,), (2) 



here Kq is a constant. This model has the advantage that it can be analytically solved |8l4l3l. ll8l420l |. 

EXACT ANALYTICAL SOLUTION 

We start with a particle moving on any of the two diabatic curves. The problem is to calculate the probability that 
the particle will still be on any one of the diabatic curves after a time t. We write the probability amplitude as 

*<-" = («M))- (^' 

where ipi{x,t) and ip2ix,t) are the probability amplitude for the two states. ^{x,t) obey the time dependent Schro 
dinger equation (we take h = 1 here and in subsequent calculations) 

9*(a;,i) ,,,, , 
i—^ ^ H'i'{x,t). (4) 

H is defined by 

fH,{x) V{x) \ 
^-[ Vix) H2{x) ) ' ^^^ 

where Hi{x) is 

F.(.)^-^^ + T..(.x). (6) 

We find it convenient to define the half Fourier Transform ^(cj) of ^{t) by 

*(w) = / *(i)e'"*dt. (7) 

Half Fourier transformation of Eq. ([4]) leads to 

l'(w) = iG(w)*(0), (8) 

where G{io) is defined by 

(w - H)G{u:) = I. (9) 

In the position representation, the above equation may be written as 



where G{x,xq]lo) is 



Writing 



^{x,uj) = i G{x,xo;uj)^ixQ,uj)dxo, (10) 



G{x,xo;u;)^{x\{w-H)-^\xo). (11) 



\G2i{x,Xo;uj) G22{x,Xq:uj) J 
and using the partitioning technique [23] we can write 

G^^ix,XQ;u) ^ {x\[u - Hi -Viu - H2)-'V]-^\xo). (13) 

The above equation is true for any general V. This expression simplify considerably if ^ is a delta function located 



-.Of .,A , K^ G^Jx, Xc; uj)G'^{xc, Xc; uj)G^{xc, xo ;uj) 

1 - K^G°i{xc, Xc;uj)G'^{xc, Xc;uj) 



Gii{x,xo;uj)^G'({x,xo;uj)-\ ^tjtWt vFJoT \ ' i^^) 
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FIG. 1: Schematic diabatic potential energy curves illustrating the model. 



where 



G°,{x,xo;Lo)^{x\iu:-H^)-'\xo), 



(15) 



and corresponds to propagation of the particle starting at Xq on the second diabatic curve, in the absence of coupling 
to the first diabatic curve. Using the same procedure one can get 



n (r,, r,, ., ,\ — KoG°{x,x^:u)G2{xa,xo;ui) 
^12[X,Xo,UJ) - i_K-^G'l(xo,x^;Lo)G'i(x,,Xa;uj} ' 



(16) 



Similarly one can derive expressions for 6*22(2^, a^o J '^) and 6*21(2;, a^oj'^)- Using these expressions for the Green's 
function in Eq. ([8]) we can calculate ^(oj) explicitly. 

The expressions that we have obtained for '^{ui) are quite general and are valid for any Vi{x) and V2(a;). However, 
their utility is limited by the fact that one must know Gi{x, xq; w) and G^{x, xq; w). It is possible to find G^(a;, xq', tj) 
only in a few limited cases and the Morse oscillator is one of them [23l425| . 



ELECTRONIC ABSORPTION SPECTRA AND RESONANCE RAMAN EXCITATION PROFILE 

CURVE CROSSING INDUCED DISSOCIATION 



In this section we apply the method to the problem involving a Harmonic oscillator and a Morse oscillator, coupled by 
Dirac delta function. We consider a system of three potential energy curves, ground electronic state and two 'crossing' 
excited electronic states (electronic transition to one of them is assumed to be dipole forbidden and while it is allowed 
to the other) i26i, 275]. We calculate the effect of 'crossing' induced dissociation on electronic absorption spectra and 
on resonance Raman excitation profile. The propagating wave functions on the excited state potential energy curves 
are given by solution of the time dependent Schrodinger equation 



d ( il}-"^^{x,t) 



'dt \r2\x,t) 



^21 (a;) Hyit^e2{x) 



ri^'ix^t) 
r2\x,t) 



(17) 



In the above equation Hmb^ei{x) and Hmb,e2{x) describes the vibrational motion of the system in the first electronic 
excited state (allowed) and second electronic excited state (forbidden) respectively 



H„ 



ib.el 



{x) 



and 



H, 



vib,e2 



(x) 



1 d^ 
1 d^ 

2m dx^ 
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(18) 



(19) 
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FIG. 2: Calculated electronic absorption spectra with coupling (solid line) and without coupling (dashed line). 



In the above ra is the oscillator's mass, l^a is the vibrational frequency of the first electronic excited state and Dp is the 
dissociation energies of the forbidden states and x is the vibrational coordinate. Shifts of the vibrational coordinate 
minimum upon excitation are given by a and &, and V12 (V21) represent coupling between the two diabatic potentials 
which is taken to be 



where Kq represent the strength of the coupling. The intensity of electronic absorption spectra is given by [27 

iGix,xo;Lo + iT)M'-,'\xo)], 
where 



G{x,xo;uj + ir) = (x|[(wo/2 + u - Ueg) + iT - Hy^bA ^\xo). 



and 



H, 



vib.e 



Hy^b^el{x) Ko\Xc){Xc\ 
Ko\Xc){Xc\ Hyib^e2{x) 



Here, F is a phenomenological damping constant which account for the life time effects. ^" (a;, 0) is given by 



^T\x,o) = 



Xt{x) 




(20) 

(21) 
(22) 

(23) 
(24) 



where Xi(x) is the ground vibrational state of the ground electronic state, ujq is the vibrational frequency on the 
ground electronic state, ea is the energy difference between the excited (allowed) and ground electronic state, and 
for the forbidden electronic state it's value is ep- Similarly resonance Raman scattering intensity can be expressed in 
terms of Green's function and is given by [27J, |28[. 



Iniio) ex I /_°1 dx /.°1 dxoWf* (x, 0) 
iG{x,xo;uj + ir)-i'r''{xo,0)\\ 



In the above VP^ (x, 0) is given by 



^f\x,0) 



Xfi^) 




(25) 



(26) 



where xAx) is the final vibrational state of the ground electronic state. As ^^{x, xq, oj) for the harmonic potential is 



23|, we can calculate G{x, xq; to). We use Eq. (|25l) to calculate the effect of curve crossing induced dissociation 



known 

on resonance Raman excitation profile. 




0.00 

8000 9000 10000 11000 12000 13000 14000 

Wavcnumbcr (cm" ) 



FIG. 3; Calculated resonance Raman excitation profile for excitation from the ground vibrational state to the first excited 
vibrational state, with coupling (solid line) and without coupling (dashed line). 



Results using the model 



In the following we give results for the effect of curve crossing on electronic absorption spectrum and resonance Raman 
excitation profile in the case where one dipole allowed electronic state crosses with a dipole forbidden electronic state 
as in Fig. [T] As in 27^, the ground state curve is taken to be a harmonic potential energy curve with its minimum at 
zero. The curve is constructed to be representative of the potential energy along a metal-ligand stretching coordinate. 
We take the mass as 35.4 amu and the vibrational wavenumber as 400 cm~^ 27| for the ground state. The first 



diabatic excited state potential energy curve is displaced by 0.1 A and is taken to have a vibrational wavenumber 
of 400 cm~^. Transition to this state is allowed. The minimum of the potential energy curve is taken to be above 
10700 cm~^ of that of the ground state curve. The second diabatic excited state potential energy curve is taken to 
be an un-displaced excited state. On that potential energy curve, the vibration is taken to have same wavenumber 
of 400 c'm~^. Its minimum is 10800 cm~^ above that of the ground state curve. Transition to this state is assumed 
to be dipole forbidden. The two diabatic curves cross at an energy of 10804.1 cm~^ with Xc = —0.02477 A. Value of 
Kq we use in our calculation is Kg — 5.54275 x 10~^^ erg. A. The lifetime of both the excited states are taken to be 
450 cm" ^. The calculated electronic absorption spectra is shown in Fig. [2l The profile shown by the dashed line is in 
the absence of any coupling to the second potential energy curve. The full line has the effect of coupling in it. The 
calculated resonance Raman excitation profile is shown in Fig. [3l The profile shown by the full line is calculated for 
the coupled potential energy curves. The profile shown by the dashed line is calculated for the uncoupled potential 
energy curves. It is seen that curve crossing effect can alter the absorption and Raman excitation profile significantly. 
However it is the Raman excitation profile that is more effected. 



CONCLUSIONS 



In our earlier paper we have proposed an exactly solvable model for the two state curve crossing problem. In this 
paper we have extended our model to deal with the case of curve crossing induced dissociation. We have analyzed 
the effect of curve crossing on electronic absorption spectrum and Resonance Raman excitation profile for a model 
consisting of a Harmonic oscillator and a Morse oscillator, coupled by Dirac delta function. We find that the Raman 
excitation profile is affected much more by the crossing than the electronic absorption spectrum. The same procedure 
is also applicable to the case where 5 is a non-local operator, and may be represented by S" = |/ > Kq < g\, f and 
g are arbitrary acceptable functions. Choosing both of them to be Gaussian will be an improvement over the delta 
function coupling model, S can also be a linear combination of such operators. 
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